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ABSTRACT 



We explore the effects of element diffusion due to gravitational settling and thermal 
and chemical diffusion on the pulsational properties of DA white dwarfs. To this end, 
we employ an updated evolutionary code coupled with a pulsational, finite difference 
code for computing the linear, non-radial g-modes in the adiabatic approximation. 
We follow the evolution of a 0.55 M0 white dwarf model in a self-consistent way 
with the evolution of chemical abundance distribution as given by time dependent 
diffusion processes. Results are compared with the standard treatment of diffusive 
equilibrium in the trace element approximation. Appreciable differences are found 
between the two employed treatments. We conclude that time dependent element 
diffusion plays an important role in determining the whole oscillation pattern and the 
temporal derivative of the periods in DAV white dwarfs. 

In addition, we discuss the plausibility of the standard description employed in 
accounting for diffusion in most of white dwarf asteroseismological studies. 
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1 INTRODUCTION 

' Asteroseismology is a method to extract information about 
, the internal structure and evolution of stars by means of 
■ the study of their oscillatory pattern. This technique, very 
sophisticated in the case of the Sun, has also undergone a 
, strong development in other stars, in particular the pulsat- 
' ing white dwarfs (WDs) (for reviews, see, e.g.. Brown & 
Gilliland 1994; Gautschy & Saio 1995; 1996). 

Pulsating WDs show multi-periodic luminosity varia- 
tions in three ranges of effective temperatures (Tcfr ) corre- 
sponding to the currently called DOV, DBV and DAV (see, 
e.g., the reviews by Winget 1988 and Kepler & Bradley 
1995). Of interest for this work are the DAVs (hydrogen- 
dominated atmospheres), or ZZ Ceti, that pulsate in the 
instability strip corresponding to 12500 K > Tofi > 10700 
K. The periodicities in the light curves of pulsating WDs 
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are naturally explained in terms of non-radial g-modes of 
low harmonic degree {I < 2), driven by the "k mechanism" 
working in a partial ionization region near the stellar surface 
(Dolez & Vauclair 1981, Winget et al. 1982). Other physi- 
cally plausible mechanism for overstability of g-modes in ZZ 
Ceti stars is the "convective driving mechanism" (see Brick- 
hill 1991 and Goldreich & Wu 1999 for details). The periods 
(P) are found within a range of 100 s < P < 1200 s and 
photometric amplitudes reach up to 0.30 magnitudes. 

Asteroseismology of WDs has recently reached impor- 
tant success, supplying independent constraints to several 
structural quantities. As a few examples we mention the 
cases of DOV PG 1159-035 (Winget et al. 1991), DBV GD 
358 (Bradley & Winget 1994), and the DAVs G117-B15A 
and R548 (Bradley 1996, 1998). 

The main observable feature in WD pulsations is the 
period pattern, which can be accurately measured. Another 
important quantity is the temporal derivative of the period 
(P), that allows to measure the cooling time-scale of WDs 
and to provide constraints on the chemical composition of 
the core. In this sense, the star DAV G117-B15A is partic- 
ularly noteworthy. Its observed periods are 215.2, 271 and 
304.4 s. For the 215.2 s mode it has been possible to find 
its temporal derivative to be P = (2.3 ± 1.4) x 10~^^s s~^ 
(Kepler et al. 2000). 

As mentioned, one of the main purposes of asteroseis- 
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mology of WDs is to disentangle the observed periodic sig- 
nals in terms of the internal structure and the evolution 
of such objects. In view of the detailed available observa- 
tions, it is very important to study the pulsational proper- 
ties of DAVs in the frame of evolutionary models as physi- 
cally sound as possible. In this regard, most of existing cal- 
culations treat the chemical profile at the hydrogen-helium 
interface (the most relevant one in the context of ZZ Ceti 
pulsations) on the basis of the equilibrium diffusion in the 
trace element approximation, hereafter EDTE approxima- 
tion, (Tassoul, Fontaine & Winget 1990, TFW; Bradley, 
Winget & Wood 1992; Bradley 1996; see also Appendix). 
In the treatment used by these authors, the chemical pro- 
file may change solely as a result of changes in the state of 
ionization in the plasma (see Appendix). Thus, if the com- 
positional transition region occurs at thermodynamical con- 
ditions at which the plasma is fully ionized, then such a 
treatment predicts fixed profiles. However, even in the case 
of thick hydrogen envelopes (when both hydrogen and he- 
lium are completely ionized deep in the star) element dif- 
fusion modifies the chemical abundance distribution within 
the star, and this is true even during evolutionary stages cor- 
responding to the ZZ Ceti domain (see Iben & McDonald 
1985, particularly their Fig. 4). 

It is the aim of this work to perform new pulsation cal- 
culations in DA WDs by relaxing both the trace element 
approximation and the diffusive equilibrium assumption. To 
this end, we carry out time-dependent diffusion calculations 
for a multicomponent plasma in a self-consistent way with 
stellar evolution. Detailed diffusion calculations consistent 
with stellar evolution have recently been performed by Mc- 
Donald, Hernanz & Jose (1998) to study the problem of 
carbon pollution in cool WDs and by Althaus, Serenelli & 
Benvenuto (2001) to assess the role played by diffusion in the 
occurrence of hydrogen thermonuclear fiashes in low-mass, 
helium-core WDs. In addition, Dehner & Kawaler (1995) 
have considered time-dependent diffusion in evolving hot 
WDs in the interest of exploring the possibility of an evo- 
lutionary link between DO PG 1159 stars and the much 
cooler DB WDs. In the context of pulsations, the change 
in the chemical composition (particularly at the hydrogen- 
helium interface) induced by diffusion processes is expected 
to affect the shape of the Ledoux term B and hence the 
Brunt- Vaisala frequency (see Brassard et al. 1991 for a dis- 
cussion of the calculation of the Brunt- Vaisala frequency in 
the context of WDs). 

In order to gauge the actual importance of time- 
dependent diffusion in the computation of theoretical P and 
P we have to calculate the WD cooling considering a full 
evolutionary code considering diffusion coupled to a pulsa- 
tional code. To our knowledge, this is the first time that 
such kind of calculation has been undertaken in the context 
of DA WDs. Note that in such kind of treatment the internal 
chemical profile is the consequence of realistic evolutionary 
models. Here, we present calculations of linear, adiabatic, 
non-radial pulsations of DAV models with a mass of 0.55 
Mq (which is representative of the mass of G117-B15A). 
In particular, we shall calculate two evolutionary evolution- 
ary sequences, one considering time-dependent element dif- 
fusion, and another one in the frame of the standard EDTE 
approximation . 

The remainder of this paper is organized as follows. In 



Section 2 we describe our evolutionary-pulsational computer 
code, paying special attention to the method for simulat- 
ing the diffusion of elements in a time-dependent approach. 
Section 3 is devoted to presenting the calculations we per- 
formed. Finally, in Section 4 we discuss our results and make 
some concluding remarks. 



2 OUR COMPUTER CODE 

2.1 Evolutionary code and diffusion equations 

The evolutionary code we employed is detailed in Althaus 
& Benvenuto (1997, 1998). This code is based on a very 
detailed and up-to-date physical description such as OPAL 
radiative (Iglesias & Rogers 1996) and molecular (Alexan- 
der & Ferguson 1994) opacities. The equation of state is an 
updated version of that of Magni & Mazzitelli (1979). High- 
density conductive opacity and neutrino emission rates are 
taken from the works of Itoh and collaborators (see Althaus 
& Benvenuto 1997 for details). Also, a complete network of 
thermonuclear reaction rates corresponding to the proton- 
proton chain and the CNO bi-cycle is included. Nuclear re- 
action rates are from Caughlan & Fowler (1988) and electron 
screening is treated as in Wallace, Woosley & Weaver (1982). 

Gravitational settling, and chemical and thermal diffu- 
sion have been fully taken into account following the treat- 
ment for multicomponent gases presented by Burgers (1969). 
Thus, we avoid the trace element approximation usually in- 
voked in most WD studies. Radiative levitation, which is 
important for determining photospheric composition of hot 
WDs (Fontaine & Michaud 1979) has been neglected. 

As a result of gravity, partial pressure, thermal gradi- 
ents and induced electric fields (we neglect stellar rotation 
and magnetic fields) the diffusion velocities in a multicom- 
ponent plasma satisfy the set of A'^ — 1 independent linear 
equations (Burgers 1969) 



dpi _ fhdp 
dr p dr 



liZicE = Kij {wj — Wi) 
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rrij Ti — rrii rj 
' rrii + rrij 
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and heat flow equation (A'' equations) 
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In these equations, Pi, pi, rii, Zi and rrii means, respec- 
tively, the partial pressure, mass density, number density, 
mean charge and mass for species i {N means the num- 
ber of ionic species plus electron). The quantities T and fcg 
are the temperature and Boltzmann constant. The unknown 
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variables are the diffusion velocities with respect to the cen- 
tre of mass, Wi, and the residual heat flows (for ions and 
electrons). In addition the electric field E has to be deter- 
mined. The resistance coefficients {Kij, Zij, zlj and z'.j) axe 
from Paquette et al (1986). 

The set of equations is completed by using the condi- 
tions for no net mass flow with respect to the center of mass 



AiUiWi = 0, 

i 

and no electrical current 
ZiUiWi — 0. 



(3) 



(4) 



In terms of the gradient in the number density we can trans- 
form Eq. (1) to 



1 



Kij (w, ~wj)+ ^ Kij Zij — 



rrii n 
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dr 



where 



d\nT 
dr ' 



(5) 



(6) 



being Ai, tuh, and g the atomic mass number, hydrogen 
atom mass and gravity, respectively. Let us write the un- 
knowns Wi, Ti and E in terms of the gradient of ion densities 
in the form 



Wi — w" — 



ions{j) 



d In rij 
dr 



(7) 



where wf* means the velocity component due to gravita- 
tional settling and thermal diffusion. With Eqs. (2) and (5) 
together with (3) and (4) we can easily find the components 
wf* and aij by matrix inversions (LU decomposition). The 
evolution of the abundance distribution throughout the star 
is found by solving the continuity equation. In particular, we 
follow the evolution of the isotopes ^H, *He, ^^C, ^"^N and 
^^O. To calculate the dependence of the structure of our WD 
models on the evolving abundances self-consistenly, the dif- 
fusion equations have been coupled to the evolutionary code. 



2.2 The pulsational code 

In order to compute the g- modes of the WD models we have 
coupled our evolutionary code to our new, finite difference, 
pulsational code described in Corsico & Benvenuto (2002), 
which solves the equations for linear, adiabatic, non-radial 
pulsations (Unno et al. 1989) . 

We describe now how these codes work together. To be- 
gin with, an interval in P and Tcff (Tcff-strip) is chosen. The 
evolutionary code computes the model cooling until the hot 
edge of the Tea-strip is reached. Then, the program calls 
the pulsation routine beginning the scan for modes. When 
a mode is found, the code generates an approximate solu- 
tion which is iteratively improved to convergence and stored. 



This procedure is repeated until the period interval is cov- 
ered. Then, the evolutionary code generates the next stellar 
model and calls pulsation routines again. Now, the previ- 
ously stored modes are taken as initial approximation to the 
modes of the new stellar model and iterated to convergence. 
Such a procedure is automatically repeated for all evolution- 
ary models inside the Tes-strip. The computational strategy 
described above has been successfully applied in fitting the 
observed periods of G117-B15A to impose constraints on 
the mass of axions (Corsico et al. 2001a) and in computing 
the period structure of low mass, helium WDs (Corsico & 
Benvenuto 2002). 

We have tested our pulsational code with two carbon- 
oxygen DA WD models of 0.5 M© and 0.85 Mq, the struc- 
ture of which was computed with the WDEC evolutionary 
code. The vibrational properties of such models were pre- 
viously analyzed by Bradley (1996). In the interests of a 
detailed comparison, we have considered a large amount of 
modes and we found that the differences between the two 
sets of modes remain below ~ 0.1%. 



3 COMPUTATIONS 

We have evolved a 0.55 Mq WD model with an internal 
carbon - oxygen chemical profile corresponding to that cal- 
culated by Salaris et al. (1997). Such a model has hydro- 
gen and helium mass fractions of Afn/Af* = 10^* and 
Muc/Mt = 10~^ respectively. These values are in good 
agreement with evolutionary predictions and are also very 
similar to those found by Bradley (1998) for the case of 
G117-B15A. The internal chemical profile of our model is 
shown in Fig. 1. It is important to mention that at the bot- 
tom of the hydrogen envelope of our model, hydrogen and 
helium are fully ionized and this is so throughout the entire 
evolutionary stages we study in the present paper. Thus, the 
chemical abundance profile predicted by the trace element 
approach remains fixed during evolution. In computing ra- 
diative opacities, we have assumed Z — 0. We have treated 
convective transport in the frame of the MLS version of the 
mixing length theory. The MLS prescription, characterized 
by a high convective efficiency, assumes the mixing length 
to be two times the local pressure scale height (see Tassoul 
et al. 1990). 

A realistic starting model for our evolutionary se- 
quences was obtained by artificially brightening an initial 
WD configuration (see Benvenuto & Althaus 1998) up to 
log L/Lq = 2. Such a procedure is known to produce an ini- 
tial sequence of some unphysical models, but then relaxes to 
the correct cooling sequence (see Althaus & Benvenuto 2000 
for further discussion) far before reaching the DAV instabil- 
ity strip. Since then on element diffusion is incorporated. 
When the model reaches Tcff = 14000 K we start pulsa- 
tional calculations. Specifically we have calculated dipolar 
{£ = 1) modes (which are usually encountered in ZZ Ceti 
light curves) with radial orders k — 1, ■ ■ ■ ,21 which cover 
a period interval 100 s < P < 1000 s. Calculations are 
stopped at Tcff — 10000 K, thus, the Tcff -strip amply em- 
braces the observed DAV instability strip. For the modes we 
have found to fulflU such conditions, we have computed peri- 
ods and eigenfunctions. For computing the Brunt - Vaisala 
frequency, we have employed the appropriate prescription 
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for degenerate models, given in Brassard et al. (1991). After 
period assessment we compute P by numerical differentia- 
tion. 



4 RESULTS AND IMPLICATIONS 

We begin by examining Fig. 2, in which we show the evo- 
lution of the chemical profile for helium resulting from time 
dependent element diffusion together with the chemical pro- 
file arising from the EDTE approximation. Note that in the 
latter case, the profile remains unchanged throughout the 
evolution because, as we mentioned, ionization is complete 
at such deep layers p|. By contrast, in the case with time- 
dependent diffusion, the chemical abundance distribution 
evolves appreciably during the ZZ Ceti evolutionary stages. 
Note that the shape of the profile in both treatments turns 
out to be markedly different particularly at the centre of 
the transition. As we shall show below, this will have an ap- 
preciable inffuence on the P and P values for some of the 
modes. 

In Fig. 3, we depict the resulting Ledoux term (panel 
A) and the squared Brunt- Vaisala frequency (panel B) at 
two selected Toh values. Let us remind the reader that the 
term B depends not only on the shape of the chemical profile 
but also on the thermal and mechanical structure of the star 
(see Brassard et al. 1991). Thus, even in the EDTE approach 
the term B changes with cooling. Because of the fact that in 
the Ledoux term there appears the derivative of the chemical 
profile (see Eq. 35 of TFW), a slight change in the slope of 
the hydrogen- helium interface translates into a noticeable 
change in B. Thus, it is not surprising the B term exhibits 
a sharp peak in the case of the trace element treatment, in 
contrast with the more physical treatment as given by non- 
equilibrium diffusion. In turn, this feature is reffected in the 
Brunt- Vaisala frequency. As it has been exhaustively shown 
by Brassard et al. (1992ab), the shape of the Brunt- Vaisala 
frequency at the chemical interfaces plays a key role in fixing 
the structure of the period pattern (e.g., mode trapping) of 
ZZ Ceti stars. 

Now let us turn our attention to the computed pulsa- 
tional modes. In Figs. 4 and 5 we show P and P corre- 
sponding to modes with I = l,fc = l,---,6 for models with 
time-dependent element diffusion and with the EDTE ap- 
proximation at the hydrogen-helium interface as a function 
of Teff . From a close inspection of these figures it can be real- 
ized that the effects of time-dependent element diffusion are 
indeed non-negligible in P and P, although for some modes 
results are very similar. We want to mention that the same 
trend has been found in modes of higher radial order (not 
shown here for brevity). Note that for the modes analyzed 
in these figures the greatest relative differences encountered 
are « 20% for P and 5% for P. 

It is worth mentioning that the differences cited above 
arise mainly from the very different shape of the interface 
profile resulting from the two treatments of diffusion inves- 
tigated here. In particular, the differences in the helium pro- 

^ In the calculations presented here, the EDTE approximation 
has been applied only to the hydrogen-helium interface, which is 
the most relevant in the context of DA WD pulsations. 



file for Xue > 0.5 (see Fig. 2) in these treatments are the 
main reason why the periods and period derivatives become 
different. In addition, there is a small contribution to the 
differences in P and P due to the evolution of the profiles 
of each chemical interface in response to non-equilibrium 
diffusion. 

Thus, as models with time dependent element diffusion 
are more physically plausible, these should be taken into 
account when an asteroseismological fit to observed periods 
is performed. Also, as P is modified, this approach should 
also be taken into account at using observed P values to 
infer the composition of the WD core. 

Before closing the paper, we would like to discuss at 
some length a major issue raised by our referee. Indeed, in 
his (her) report, our referee asked us to look for the underly- 
ing physical reasons for the differences between the standard 
treatment of EDTE approximation and our full treatment 
of time dependent element diffusion. To be specific, the ref- 
eree asked us whether such differences are due mostly to the 
relaxation of the trace element approximation or the equi- 
librium hypothesis. 

In order to find the answer, we decided to perform a 
simple immerical experiment: we simulate the equilibrium 
diffusion conditions with our full code. Equilibrium diffusion 
would be a good approximation if the diffusion timescale 
were much shorter than the evolutionary one. Then, in order 
to simulate this situation, we simply assumed the diffusion 
time step to be several times the evolutionary one^ This is 
equivalent to assuming that the whole diffusive process oc- 
curs several times faster. We computed the evolution of our 
WD model under this hypothesis. The result was that the 
WD did not evolve along the cooling branch hut instead suf- 
fered from a hydrogen thermonuclear flash. Physically, the 
reason for this behaviour is that, if diffusion were plenty of 
time to evolve to equilibrium profiles, then the tail of the 
hydrogen profile would be able to reach hot enough layers 
to be ignited in a fiash fashion. The fact that at imposing 
equilibrium diffusion conditions the star undergoes a hydro- 
gen thermonuclear flash, while in the detailed, self-consistent 
time dependent diffusive treatment the star cools down qui- 
escently, clearly shows the incorrectness of the hypothesis 
of equilibrium. In view of this, we are forced to conclude 
that there is no way other than to abandon the idea that 
the shape of the internal proflles in the WD is determined 
by equilibrium diffusion. This conclusion is valid at least for 
massive hydrogen envelopes. We think that the only physi- 
cally sound way to compute such proflles, a key ingredient 
in asteroseismological studies, is to calculate the WD evo- 
lution in a self-consistent way with time dependent element 
diffusion and nuclear burning. 

In addition, some words are in order about the standard 
treatment of equilibrium diffusion, which is based on the 
work of Arcoragi & Fontaine (1980). Recent asteroseismo- 
logical studies of DAV WDs (Clemens 1994, Bradley 1998, 
2001) seem to favour large values for the thickness of the hy- 
drogen envelope (Mh/M« ~ 10~*). At the thermodynamical 
conditions relevant to DAV WD models, we flnd that most 
of the hydrogen-helium interface occurs at degenerate con- 



I Here, because of numerical reasons, we assumed a diffusive time 
step a hundred times the evolutionary one. 
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ditions. Because the Arcoragi & Fontaine (1980) equations 
are valid for non- degenerate conditions, we should remark 
that this treatment cannot be applied to the modeling of 
the hydrogen-helium interfaces in DA WDs with massive 
hydrogen envelopes. 

We should also remark that, in the case of the trace 
element approximation, we have found that the object does 
not undergo any thermonuclear flash. This is another arti- 
fact of the approximation, due to the ad hoc truncation of 
the profile at some low density (see Appendix). Indeed, in 
our numerical experiments, we have found that the stellar 
model evolves along the WD cooling track if we truncate 
the hydrogen profile at Xh = 10~^. However, if we allow 
the equilibrium hydrogen profile to extend to slightly lower 
abundances (e.g., Xh ~ 10"'') the model experiences a ther- 
monuclear flash! 

The results presented in this paper indicate that a more 
extensive and systematic exploration of asteroseismology of 
DAV in the frame of detailed evolutionary models consid- 
ering time dependent element diffusion is worth being done 
and it will be the subject of further papers. While the present 
paper was in process of reviewing, some interesting results 
about the effects of diffusion on mode trapping have been 
presented in Corsico et al. (2001b). 
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For the sake of completeness, we describe here the equations 
employed to include diffusion processes in most of WD pul- 
sational studies. This approximation is based on the work 
of Arcoragi & Fontaine (1980) (see also Tassoul et al. 1990). 
Here we limit ourselves to comment on the most important 
aspects involved in this treatment. 

To begin with, Arcoragi & Fontaine (1980) assume a 
stellar plasma made up of two-ionic species with average 
charge Zi and Z2 and atomic weight Ai and A2. In addi- 
tion, thermal diffusion is neglected and an ideal gas equation 
of state is considered under the assumption that the plasma 
is sufficiently diluted. Under these approximations, the dif- 
fusion velocity u!i2 reads 



Wx2 



dr 

A2Z1 -A1Z2 eE ' 
Ai+'yA2 ksT 



dlnc2 ^ A2 — A\ d\np ^ 



Ar + 'yA2 dr 



(8) 



D12 is the diffusion coefficient, and d, the number concen- 
tration of ions of species i, is defined by 

_ rii Pi 



ni +n2 pi + P2 



(9) 
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being pi the partial pressure. E is the electric field, given by 



''^""'''^^l(^l + l)+Z2(Z2 + l)7' 

and 7 is defined as 

^2 _ P2 _ £2 

ni pi ci ' 



7 



(10) 



(11) 



The remainder of the symbols have the usual meaning. 
Notice that Eqs. and ( |ll| ) are valid in an isothermal 
medium, i.e., we are neglecting temperature gradients. 

Now, we impose equilibrium diffusion, by assuming 
Wi2 = 0, and from Eqs. - we get the ordinary dif- 
ferential equation for the equilibrium profile (Eq. A5 of Ar- 
coragi & Fontaine 1980). In the trace element approximation 
(7 ^ 1), for the specie 2 considered as a trace, we get 



a In 



C2 



dr 
where 



a2- 



dlnq 
dr ' 



A. 



(12) 



(13) 



02 = -r^ (l + Zi) -Z2-I. 

For the purpose of application, Tassoul et al. (1990) 
divide the transition zone into two parts: an upper one in 
which element 1 is dominant and element 2 is a trace, and 
a lower one in which the roles of the respective elements 
is reversed. For the upper region the abundance profile of 
element 2 considered as a trace is given by Eq. (|l^) and 
for the lower part of the transition zone the abundance of 
element 1 (considered as a trace) is given by 

9 In ci 9 In g 

= ai- 



dr 
where 

Ql 



dr 



(14) 



(15) 



q is the mass fraction (1 — Mr/M,). The integration of Eqs. 
( p^ and ( p^ ) gives the equilibrium abundance profiles: 



C2 



29 



and 

i<3 



(upper region of interface) 



(lower region of interface) 



(16) 



(17) 



By invoking the condition of continuity in the middle 
point of the inteface, we obtain the relation 

feg^^ = fci9m = (18) 

where is the mass fraction where the abundances of two 
species are equal; the qm value is obtained by forcing the 
mass conservation of element 1. Thus, in the case of the 
hydrogen- helium transition region, the outer mass fraction 
of hydrogen {qn = M-n/M^) is employed for computing qm- 
Note that possible changes in the equilibrium profiles re- 
sult only from slight changes in the ionization states of the 
elements present at the interfaces, i.e., variations in the ex- 
ponents Ql and 0:2 (Tassoul et al. 1990). 

To implement this approach to modeling of the 
hydrogen- helium transition zone, it is necessary to set small 
abundances to zero in order to avoid having a tail of hydro- 
gen in regions deep enough where carbon is abundant (for 
this case we should generalize the above treatment for three 



species). Moreover, if we do not do so, hydrogen would be 
present at layers hot enough to force the star to undergo a 
thermonuclear fiash. This is the case, at least for thick hy- 
drogen envelopes like those we have treated here which, in 
turn, are the favoured ones by current asteroseismological 
studies (see the main text). 
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Figure 1. The internal chemical profiles of the 0.55 Mq carbon-oxygen WD model for hydrogen, helium, carbon and oxygen at an 
effective temperature of 14000K. In the case of the EDTE approximation, the fixed profiles are represented by dotted lines. Profiles for 
models in which time dependent element diffusion has been considered are represented by solid lines, q is the outer mass fraction defined 
by g = 1 - Mr/M,. 
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Figure 2. Panel A: The internal chemical profile for helium in terms of the outer mass fraction q at the hydrogen-helium interface. Solid 
line corresponds to a model at a temperature of 14000K in which time-dependent diffusion is considered. Long dashed line means the same 
treatment but for a model at lOOOOK. Finally dotted line corresponds to the EDTE approximation prediction. Panel B: The chemical 
profile for helium as given by the EDTE approximation according to a 0.50 MQmodel calculated by Bradley (private communication). 
Note that the shape of the profile in both models with diffusive equilibrium is the same. 
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Figure 3. In the upper panel it is shown the Ledoux term B at the hydrogen-helium interface in the case of time-dependent diffusion at 
the Tgfj values of 14000K and lOOOOK as given by solid and long dashed lines, respectively. Dotted and dot-dashed lines correspond to 
the same temperature values but for the EDTE approximation. In the lower panel the square of the Brunt- Vaisala frequency is shown 
for the same cases analyzed in panel A. For details, see text. 
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Figure 4. Period and period derivative for £ = 1, fc = 1, 2, 3 modes for a 0.55 Mq carbon-oxygen WD model in at a T^fj interval containing 
the DAV instability strip. Solid lines (filled squares) correspond to periods (period derivatives) computed considering non-equilibrium 
diffusion while dotted lines (empty squares) depict periods (period derivatives) computed according to the EDTE approximation at the 
hydrogen-helium interface. For discussion of the results, see text. 
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Figure 5. Same as figure 4 but for modes with I = 1, A; = 4, 5, 6. 
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